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We study the interaction between low-lying transverse col- 
lective oscillations and thermal excitations of an elongated 
Bose-Einstein condensate by means of perturbation theory. 
We consider a cylindrically trapped condensate and calculate 
the transverse elementary excitations at zero temperature by 
solving the linearized Gross-Pitaevskii equations in two di- 
mensions. We use them to calculate the matrix elements 
between thermal excited states coupled with the quasi-2D 
collective modes. The Landau damping of transverse collec- 
tive modes is investigated as a function of temperature. At 
low temperatures, the damping rate due to the Landau de- 
cay mechanism is in agreement with the experimental data 
for the decay of the transverse quadrupole mode, but it is too 
small to explain the slow experimental decay of the transverse 
breathing mode. The reason for this discrepancy is discussed. 

03.75.Fi, 02.70.Lq, 67.40.Db 



I. INTRODUCTION 

Transverse collective oscillations have been recently ex- 
cited in the radial plane of an elongated cylindrically 
symmetric condensate They are quasi-2D excita- 

tions and indeed, at very low temperatures, the measured 
frequencies of the transverse monopole and quadrupole 
modes are found to be in very good agreement with the 
theoretical predictions obtained by linearizing the Gross- 
Pitaevskii equation for a trapped two-dimensional (2D) 
system. Moreover, it has been experimentally found that 
the transverse breathing mode of an elongated conden- 
sate exhibits unique features. First, the excitation fre- 
quency is close to twice the radial trapping frequency, it 
is almost independent of the strength of the two-body 
interaction and the number of particles, and it is nearly 
independent of temperature. And second, the transverse 
breathing mode has a very small damping rate compared 
to other modes. The measured damping rate of the trans- 
verse quadrupole mode is approximately one order of 
magnitude larger than the one of the transverse breath- 
ing mode. This is in contrast to the 3D case, where 
both monopole and quadrupole modes have similar decay 
rates. These striking properties experimentally found for 



the transverse breathing mode are in agreement with the 
unique features predicted for the breathing mode in 2D 
isotropic condensates [^|-§}: its frequency is a universal 
quantity and it is an undamped mode. One can conclude 
that these features are characteristic of the reduced di- 
mensionality thus, providing a demonstration of the two- 
dimensional nature of the transverse oscillations of an 
elongated condensate. Therefore, a highly anisotropic 
cigar-shaped trap can be considered as a first approxi- 
mation as an infinite cylinder along the ^-direction. 

In the present paper we consider a cylindrical con- 
densate and calculate at zero temperature the trans- 
verse spectrum of excitations, by solving the linearized 
Gross-Pitaevskii equation in 2D, and the modes with 
non-zero momentum k along the longitudinal axis, as 
well. The main purpose is to investigate the coupling 
between quasi-2D (transverse) collective oscillations and 
thermal excitations in the collisionless regime. We study 
the Landau damping as a decay mechanism of transverse 
collective modes as a function of temperature, in which 
a transverse collective oscillation is annihilated in a col- 
lision with a thermal excitation to give rise to another 
excitation. 

This paper is organized as follows. In Sec. II we calcu- 
late the transverse spectrum of excitations of a cylindri- 
cal condensate by solving the linearized Gross-Pitaevskii 
equation in 2D within the Bogoliubov theory. We obtain 
also the modes with non-zero k. In Sec. Ill we briefly re- 
call the perturbation theory for the interaction between 
collective modes of a condensate and thermal excitations 
(developed in Ref . || ) . We apply this theory to calculate 
the decay rate due to Landau damping of a transverse os- 
cillation. In Sec. IV we discuss the main results for the 
transverse quadrupole mode. The situation with the sin- 
gular transverse breathing mode is also discussed. 

II. ELEMENTARY EXCITATIONS OF A 
CYLINDRICALLY TRAPPED CONDENSATE 

At low temperature, the elementary excitations of a 
trapped weakly interacting degenerate Bose gas are de- 
scribed by the time-dependent Gross-Pitaevskii (GP) 



1 



equation for the order parameter: 
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where g = 4irh 2 a/M is the interaction coupling constant, 
fixed by the s-wave scattering length a, M is the atomic 
mass and the order parameter is normalized to the num- 
ber of atoms in the condensate N. The confining poten- 
tial is usually cylindrically symmetric and it is given by 
V ext (r x ,z) = M(u\r\ + uj 2 z z 2 )/2, with r\ = x 2 + y 2 . 
The ratio between axial (u> z ) and radial trapping 
frequencies defines the anisotropy parameter of the trap 
A = lj z /lj±. 

Experimentally, in Ref. jlj , the transverse modes have 
been excited in the radial plane of a highly anisotropic 
cigar-shaped trap, with A ~ 0.0646. Thus, as a reason- 
able approximation, we consider an idealized cylindrical 
trap that is uniform in the z-direction (lo z — 0) and has 
isotropic trapping potential in the radial plane 
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The harmonic trap frequency lo±_ provides a typical 
length scale for the system, a± — {fi/Muj^) 1 / 2 . We as- 
sume a cylindrical condensate with large longitudinal size 
L compared to the radial size R, and with a number of 
atoms per unit length N± = N/L. 

At low temperatures, the dynamics of the condensate 
is described by the linearized time-dependent GP equa- 
tion. Assuming small oscillations of the order parameter 
around the ground-state value 



*(r,t) = e 
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where fi is the chemical potential, the normal modes of 
the condensate can be found by seeking fluctuations such 
as 



5V(r,t) =u(r)e- luJt + v*(r)e i 



(4) 



The functions u and v arc the "particle" and "hole" com- 
ponents characterizing the Bogoliubov transformations, 
and w is the frequency of the corresponding excited state. 

The transverse excitations of a highly elongated con- 
densate can be obtained by neglecting oscillations along 
the longitudinal direction and thus, looking for fluctua- 
tions in the transverse plane of the form 



(5) 



6V(r±,6,t) = U{r ± )e" n6 e- luJt + V*(r ± )e 



where we have used cylindrical coordinates (r_i_, 0,z). 

Transverse excitations are quasi-2D modes. To calcu- 
late Landau damping we must also consider thermal ex- 
citations having a quasi-continuous spectrum at L 3> R 
with the finite longitudinal wave vector k along z direc- 
tion. To study these modes we look for solutions that 
correspond to propagating waves of the form |||| 



5*(r,i) = U k (r ± )e^ me+kz ^e- tuJt + 
Vk(r ± )e- l( - m9+kz) e wt . 
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Inserting (||) , (Q) and ^ in Eq. ([!]) and retaining terms 
up to first order in u and v, it follows three equations. 
The first one is the nonlinear equation for the order pa- 
rameter of the ground state, 



[H + g*ftr±j\ *o(r±) = n*o(r±) , 



(7) 



where H = ~{h 2 /2M)W 2 r± + V cxt (r ± ), and W 2 r± = 
d 2 /dr\ + rj 1 d/dr±; while Uk(r±) and 14(rj_) obey the 
following coupled equations M : 
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We solve the linearized GP equations Eqs. (|7[]9|) at zero 
temperature, to obtain the ground state wave function 
"to, the eigenfrequencies uo{k) and hence the energies 
E = fiLO of the excitations, as well as the correspond- 
ing functions (u, v). 

In a cylindrical trap the modes of the condensate are 
labeled by i = (n, m, k), where n is the number of nodes 
in the radial solution, m is the z-projection of the or- 
bital angular momentum, and k is the longitudinal wave 
vector. The quantum numbers n and m are discrete 
(n = 0,1,2,..., and m — 0, ±1, ±2, ....) whereas k for 
highly elongated condensates (L ^> R) is a real number, 
and thus, it is a continuous index. Moreover, it can be 
seen from Eqs. (g) and ([)]), that levels with m = ± | m | 
or k = ± | k | are exactly degenerated. 

The excitation spectrum of a cylindrical condensate 
presents different branches labeled by (m, n) that are 
continuous and increasing functions of k || . For a given 
mode, i.e., fixed n and m, the frequency of excitation is 

We have numerically checked from Eqs. (g) and (g) 
that the state with (m = 0,n = 1) and k = corre- 
sponds to the transverse breathing mode with excitation 
frequency equal to 2wj_ , and that it is independent of the 
strength of the two-body interaction potential and the 
number of particles. Of course, according to the univer- 
sal nature of this mode noted above, this general result 
is also valid in the Thomas-Fermi (TF) limit, where the 
dispersion relation can be found analytically How- 
ever, the functions u and v for this mode do not coincide 
in general with their Thomas- Fermi limit. 

The state with (m = 2,n = 0) and k = is a trans- 
verse quadrupole mode whose frequency depends on the 
two-body interaction, ranging from the TF value \[2uj^ , 
valid for large condensates (aN± > 1 or equivalently 
fi 3> hcuj_), to the noninteracting value 2uj±. Therefore, 
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the transverse quadrupole has to be obtained numerically 
by solving Eqs. @) and (J)|) in each particular value of 
the dimensionless parameter aN±. In the present work, 
we use the numerical solutions to calculate the Landau 
damping of collective modes. 

It is well known that the dipole mode (m = 1, n = 0) 
is unaffected by two-body interactions due to the trans- 
lational invariance of the interatomic force which cannot 
affect the motion of the center of mass ]j| and its excita- 
tion frequency is lu>±. It is worth stressing that the in- 
dependence of the interaction of the transverse monopole 
mode, is a unique property of 2D systems, related to the 
presence of a hidden symmetry of the problem described 
by the two-dimensional Lorentz group SO(2,l) B. 



III. LANDAU DAMPING 

Let us consider a collective mode with frequency Q osc , 
and the excited states i,j available by thermal activa- 
tion, with energies Ej and Ej , respectively. Suppose that 
this collective mode has been excited and, therefore, the 
condensate oscillates with the corresponding frequency 
Sl osc . Due to interaction effects, the thermal cloud of ex- 
citations can either absorb or emit quanta of this mode 
producing a damping of the collective oscillation. We 
want to study the decay process in which a quantum of 
oscillation hfl osc is annihilated (created) and the i-th ex- 
citation is transformed into the j-th one (or viceversa). 
The energy is conserved during the transition process, 
therefore 



Ej — Ej -I- h£l n 



(10) 



where we assume that Ej > Ej. This mechanism is 
known as Landau damping. Another possible decay 
mechanism, also due to the coupling between collective 
and thermal excitations, is the Beliaev damping [fl0|| , 
which is based on the decay of an elementary excitation 
into a pair of excitations. 

Let us define the dissipation rate 7 through the follow- 
ing relation between the energy of the system E and its 
dissipation E: 



E = -27J5. 



(11) 



Assuming that the damping is small, one can use pertur- 
bation theory to calculate the probabilities for the tran- 
sition between the i-th excitation and the j-th one, both 
available by thermal activation, yielding the following ex- 
pression for the Landau damping rate f|,|ll| : 



a, 



ij 



lij — ^osc) , 



(12) 



where the transition frequencies u>ij — (Ej — Ej)jh are 
positive, and the delta function ensures the energy con- 
servation during the transition process. We have assumed 
that the thermal cloud is at thermodynamic equilibrium 



and the states i, j are thermally occupied with the usual 
Bose factor /j = [exp(£ , j/fcsT) — 1] — 1 - The "damping 
strength" has the dimensions of a frequency and is given 
by 



lij = 717^ — I A v I (fi ~ fj)- 
A j 'osc 



(13) 



The matrix element that couples the low-energy col- 
lective mode (u osc ,v osc ) with the higher energy single- 
particle excitations (for which we use the indices is 
1 

Aij = 2g J dr ^ [(u*Vi + v*Vi + u*Ui)u osc 

+ (v*Ui + v*Vi + u*Ui)v osc }. (14) 

In this work we calculate the quantities 7y by using the 
numerical solutions u and v of Eqs. (0-^) into the inte- 
grals (fL4|), avoiding the use of further approximations in 
the spectrum of excitations. In cylindrically symmetric 
traps, i — (n, m, k) and from (^|) the longitudinal excita- 
tions are u nm k{r) = U nm k(r±) e lm0 e tkz anc j the energies 
Enmk are not degenerate. 

We are interested in the decay of transverse collective 
excitations (J|) u osc (r±,8) = U osc (rj_)e m ° acd , due to the 
coupling with thermally excited states The func- 

tions it osc and v QSC do not depend on z, and from Eq. jl^ ) 
it is straightforward to see that the matrix element 
couples only energy levels with the same quantum 
number k = fc, = kj, and Am =| rrii — rrij |= m osc . To 
calculate the damping rate 7 of the collective oscillation 
we have to sum up over all couples of excited states (i, j), 
or equivalently, over all pairs of sets (rii rrii k , nj rrij fc), 
that satisfy the above conditions. Since k is a continuous 
quantum number, the sum over k becomes and integral 
and Eq. dl|) yields 
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where the sum is restricted to pairs that satisfy Am = 
m osc . The frequency of a mode is a continuous function 
of k, therefore the transition frequency is = ujij(k). 
Then, in Eq. (15) the integration over k yields 
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where k is the wave vector in which the conservation of 
energy is verified, i.e., | u>i(k) — ojj(k) |= f2 osc , and we 
have used the known property of delta functions 



J-00 ^ I Q'(k) I 



(17) 
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IV. RESULTS 

We want to calculate the characteristic decay rates due 
to Landau damping of the transverse low-lying collec- 
tive modes of an elongated Bose-Einstein condensate as 
a function of T. We consider the collective excitations in 
the collisionless regime, which is achieved at low enough 
temperatures. 

For a fixed number of trapped atoms, the number of 
atoms in the condensate depends on temperature. At 
zero temperature the quantum depiction is negligible and 
all the atoms can be assumed to be in the condensate 
[jjjl . At finite temperature the condensate atoms coexist 
with the thermal bath. However, at low enough temper- 
atures the excitation spectrum can be safely calculated 
by neglecting the coupling between the condensate and 
thermal atoms It means that the excitation energies 
at a given T can be obtained within Bogoliubov theory 
at T = normalizing the number of condensate atoms to 
the corresponding condensate fraction at that tempera- 
ture and assuming that they are thermally occupied with 
the Bose factor. 

In order to present numerical results we choose a gas 
of 87 Rb atoms (scattering length a = 5.82 x 1CP 7 cm) 
confined in an elongated cylindrical trap with frequen- 
cies uj z w and oj± = 219 x 2ir Hz, that corresponds to 
an oscillator length a± = 0.729 x 1CP 4 cm. The num- 
ber of condensate atoms per unit length is taken to be 
Nj_ = N/L = 2800 aj 1 . These conditions are close to 
the experimental ones of Ref . H . 

We have solved the linearized GP equations to obtain 
an exact description of the ground state ^o(r±) and the 
normal modes of the condensate within Bogoliubov the- 
ory without using the Thomas- Fermi or Hartree-Fock ap- 
proximations. We have obtained the following numerical 
results: /i = 9.526 fkvj_, and the excitation frequencies 
of the transverse monopole Qm — 2.0 u>±, and of the 
transverse quadrupole mode £Iq — 1.436 hu>±. We have 
calculated the branches of the excitation spectrum of the 
cylindrical condensate, labeled by (m, n) as a function of 
k. 

To calculate the damping rate of a collective mode f2 OS c 
we have to obtain, first, all pairs of levels of the exci- 
tation spectrum that satisfy both the energy conservation 
of the transition process Ej — Ei = hfl osc ^ and the selec- 
tion rules given by the particular collective mode under 
study. And then, calculate the corresponding matrix ele- 
ments. We have restricted our calculations to levels with 
energy E < 7/i. The contribution of higher excited levels 
can be neglected since their occupation becomes negligi- 
ble in the range of temperatures we have considered. 

It is worth recalling that the dipole mode is undamped 
H since it is not affected by two-body interactions. 
We have checked that within our formalism the Landau 
damping of the transverse dipole mode is zero, as ex- 
pected. 

First of all, we have calculated the Landau damping 



of the transverse quadrupole mode (k = 0, m = 2 and 
n = 0) as a function of T, jq(T). For such a mode, 
the functions m osc and v osc depend on the radial coor- 
dinate rj_ and have also angular dependence (~ e ±l2e ). 
Therefore, from Eq. ( p^ ) it follows the selection rules 
corresponding to transverse quadrupole- like transitions: 
ArTi =| rrij — mi |= 2 and Afc = 0. Only the energy levels 
that satisfy these conditions can be coupled by the 
interaction. The integration of the radial part has to be 
done numerically. 

Fixed the number of condensate atoms, we have calcu- 
lated the Landau damping for the transverse quadrupole 
mode as a function of temperature. At zero temperature 
we have obtained the excitation spectrum within Bogoli- 
ubov theory. We have found that there are 34 pairs of lev- 
els that verify the transverse quadrupole selection rules 
and the energy condition ( |Io| ) at finite wave vectors, i.e., 
at k 7^ 0. We have calculated the matrix element Aij 
for each allowed transition, and then, the corresponding 
damping strength 7^ at a given T. Summing up the 
damping strengths for all transitions, we have obtained 
the damping rate for the transverse quadrupole oscilla- 
tion, due to the coupling with thermal excitations. 

In Figure |l| we plot the damping rate versus ksT/fi 
for the transverse quadrupole mode for a fixed number 
of condensate atoms per unit length in a cylindrical trap. 
As expected, Landau damping increases with tempera- 
ture since the number of excitations available at thermal 
equilibrium is larger when T increases. One can distin- 
guish two different regimes, one at very low T (ksT <C //) 
and the other at higher T, where damping is linear with 
temperature. However, this linear behaviour appears at 
relatively small temperature (knT ~ /x) in comparison 
to the homogeneous system || where this regime oc- 
curs at ksT ^> /x. For the highly elongated conden- 
sate of Ref. (!]], it has been measured at ksT ~ 0.7/i a 
quality factor |13| ] for the transverse quadrupole mode 
Q = Qq/(2jq) ~ 200 which gives a damping rate 
7q XP ~ 10~ 2 u>±. This value is in agreement with our 
calculations of the Landau damping for this mode (see 
Fig. 1). A more quantitative comparison is not possi- 
ble at present due to the absence of detailed measure- 
ments of the damping of this mode. However, it is rea- 
sonable to conclude that the transverse quadrupole mode 
in elongated condensates decays via Landau damping. It 
is worth noting that the order of magnitude of the damp- 
ing rate of the transverse quadrupole mode is the same as 
the one previously estimated for a uniform gas fi pT] ltjj, 
for spherical traps jllj and for anisotropic traps J17 18 
as well. 

The Landau damping for other modes can be calcu- 
lated analogously. However, the physical situation for 
the decay of the singular transverse monopole mode turns 
out to be much more complicated. Formally the calcula- 
tions can be committed in the same way. The quantum 
numbers for this mode are k = 0, m = and n = 1. The 
functions u sc arid v osc depend only on the radial coordi- 
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nate r± , and the matrix element Ay couples only energy 
levels with the same quantum numbers m and k. 

That is, the selection rules corresponding to transverse 
monopole-like transitions are Ak = (as we have already 
discussed in the previous section) and Am = 0. 

From the spectrum of elementary excitations, we have 
found that there are only 16 pairs of levels that verify the 
transverse monopole selection rules and the energy con- 
dition ( |l0| ) at k 7^ 0. The damping strengths associated 
to the allowed transitions are smaller than the jij val- 
ues associated to quadrupole-like transitions. As a result 
from our calculations the damping occurs quite small. At 
T ~ 1.5/i we obtain a Landau damping ~ 7 x 10~ 5 cj^ 
that is one order of magnitude smaller than the exper- 
imental decay jl3| 7^J P ~ 6 x 10~ 4 cl>^ measured in jj]. 
Thus, Landau damping cannot explain the experimen- 
tally observed decay of the transverse breathing mode. 

Actually, there are reasons to believe that even this 
small value of Landau damping obtained in the present 
calculations, could be larger than the true one. The point 
is that in our calculations, based on perturbation the- 
ory, the thermal cloud is assumed to be in the state of 
thermodynamic equilibrium and thus one can neglect its 
motion. This is correct for the case of quadrupole os- 
cillations, where the frequencies of the collective motion 
of the condensate and the cloud are different. However, 
due to the unique peculiar nature of the 2D breathing 
mode, the breathing oscillation of the thermal cloud has 
the same frequency 2to± as the condensate. Then, if the 
oscillation is excited by the deformation of the trap, as it 
takes place in the present experiment jlj], the cloud will 
oscillate in phase with the condensate what will further 
decrease the Landau damping [ jl9| . Therefore, this quasi- 
2D monopole mode in highly elongated condensates may 
decay via another damping mechanism, different from the 
Landau damping. 

After finishing this work it has appeared a preprint by 
Jackson and Zaremba |2^| in which the effect of the colli- 
sions between elementary excitations has been taken also 
into account. This effect is usually small in comparison 
with the Landau damping. Nevertheless, for this peculiar 
mode the collisions seems to be the dominant damping 
mechanism, yielding a decay rate that it is in agreement 
with the experimentally measured value jl]]. Note, how- 
ever, that the experiments are produced at quite high 
amplitude of oscillations and non-linear effects can play 
an important role. For example, as it was noted in [p2[ , 
the transverse breathing mode exhibits a "parametric" 
instability due to the decay into two or more excitations 
with non-zero momentum along the axis. In Ref. |23|] it 
has been considered the process in which the transverse 
breathing oscillation decays into two excitations propa- 
gating along z with momentum k and —k. This damping 
mechanism is active also at zero temperature. But in this 
case the damping rate depends on the amplitude of the 
breathing oscillation. 



V. SUMMARY 

We have investigated the decay of low-lying transverse 
oscillations of a large cylindrical condensate. First of all, 
we have calculated the transverse normal modes, as well 
as the excitations with finite k, of the condensate by solv- 
ing the linearized time-dependent Gross-Pitaevskii equa- 
tion. And then, within the formalism of Ref. H], we 
have calculated numerically the matrix elements associ- 
ated with the transitions between excited states allowed 
by the selection rules of the transverse collective modes. 
Within a first-order perturbation theory and assuming 
the thermal cloud to be in thermodynamic equilibrium, 
we have studied the Landau damping of transverse col- 
lective modes due to the coupling with thermal excited 
levels as a function of temperature. For the damping 
rate of the transverse quadrupole mode, we have found a 
good agreement with the experimentally measured value 
of Ref. (lj . Whereas our result for the transverse breath- 
ing mode is one order of magnitude smaller than the ex- 
perimental decay. One can conclude that in a highly 
elongated cylindrically symmetric condensate, transverse 
quadrupole oscillations decay via Landau damping mech- 
anism, but the transverse breathing mode which has an 
anomalously small measured damping rate, may decay 
via other damping mechanisms. 
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FIG. 1. Dimensionless damping rate of the transverse 
quadrupole mode ^q/uj± as a function of ksT/jj, for a cylin- 
drical condensate with N± — 2800a^ atoms of 87 Rb per unit 
length. 
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